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Abstract 

We examine hybrid thermal-nonthermal synchrotron spectra from a sphericaUy symmetric, optically- 
thin wind, taking into account the relativistic effect. In the relativistic flow from the central object, due to 
the relativistic beaming effect, the observed spectra often shift towards high frequency and high intensity 
directions. In the optically thin outflows, however, we find that the intensity of the observed spectra 
decreases compared with that of the emitted ones, although the peak frequency shifts towards the high 
frequency direction. This is because in the optically thin outflows we can see the far side flows that go 
away from the observer. We thus carefully consider optically thin relativistic flows around a black hole 
such as Sgr A*. 
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1. Introduction 

It is now well established that black holes with accre- 
tion disks and jets are the central engine of various active 
phenomena, such as active galaxies and quasars, X-ray bi- 
naries and micro-quasars, and gamma-ray bursts, and so 
on (cf. Kato et al. 2008). However, the nature of the 
central engine is not yet well understood, due partly to 
poor observational resolutions, and due partly to so many 
physical complications. In order to clarify the nature of 
the central engine, many researchers have examined emis- 
sion spectra and observational appearance of the relativis- 
tic flows around black holes. Depending on the mass- 
accretion and outflow rates, the relativistic flows would 
be optically thick in some cases while optically thin in 
other cases. 

In any cases, the gravitational energy released in the 
mass accretion process is converted into the tremendous 
radiation, magnetic, and kinetic energies. As a result, 
the relativistic jets and winds often blow off from these 
systems — a black hole wind. 

In the optically thick cases, Abramowicz et al. (1991) 
have pointed out that the optical depth should be care- 
fully considered in the relativistic flows. They found that 
the shape of the apparent photosphere of a spherical wind 
becomes concave at a high speed regime due to the rela- 
tivistic effect on the optical depth. 

Recently, furthermore, Sumitomo et al. (2008) has 
firstly examined the observational appearance of spheri- 
cally symmetric, relativistic massive winds (see also Fukue 
& Sumitomo 2009). They also calculated the comoving 
and observed luminosities, and found that the luminosities 
increases with the velocity but decreases with the mass- 
outflow rate. 

In the optically thin cases, however, observed spectra of 



black-hole winds has not been well examined up to now. 
in particular the relativistic effect such as Dopplcr and 
aberration effects. 

In the optically thin relativistic flows, the synchrotron 
emission is believed to be a key ingredient, and many 
researchers has discussed its role (Colpi et al. 1984; 
Begelman & Chiueh 1988; Bisnovatyi-Kogan & Lovelace 
1997; Quataert & Gruzinov 1999; Quataet & Narayan 
1999; Gruzinov & Quataert 1999; Medvedcv 2000). There 
may exist thermal and non-thermal electrons, and in or- 
der to examine the synchrotron emission we should both 
type of electrons (Zdziarski et al. 1993; Ghisellini et al. 
1993; Li et al. 1996; Ghisellini et al. 1998; Ozcl et al. 
2000). 

Of these, Ozel et al. (2000) obtained the hybrid 
thcrmal-nonthcrmal synchrotron emission spectra from 
hot accretion flows, and reproduced the observed spectra 
from the Galactic center source, Sgr A* (see also Coppi 
1999). In their calculations, they considered the hot ac- 
cretion flow, so called ADAF solutions (Narayan & Yi 
1995; Narayan et al. 1996). In addition, they did not con- 
sidered the relativistic effect in the high velocity regime, 
where the flow speed is on the order of the speed of light. 

Thus, in the present study we consider an optically thin 
hot wind from a black hole, and examine emitted and ob- 
served spectra of hybrid thermal-nonthermal synchrotron 
emissions, taking into account the relativistic effect in the 
high velocity regime. 

In section 2, the present wind model and the calculation 
method are briefly described. In section 3 the results are 
shown. Final section is devoted to concluding remarks. 
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2. Wind Model and Synchrotron Emissivity 

In this section we describe the present wind model, and 
summarize the synchrotron emissivity model. 

2.1. Wind Model 

We assume that a spherically symmetric, relativistic 
wind blows off from a central object. As a central object, 
we assume a non-rotating black hole (Schwarzschild black 
hole); the Schwarzschild radius is defined by Tg = 2GM/c^, 
where G, M, and c represent the gravitational constant, a 
black hole mass, and the speed of light, respectively. We 
use the spherical coordinates {R,9,Lp) and the cylindrical 
coordinate (r, (/3,z), whose z-axis is along the linc-of-sight 
direction. 

From continuity equation for the spherically symmetric 
stationary wind, the rest mass density po measured in the 
comoving frame varies as 
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where M is the mass-outflow rate, /3 {— v/c) the velocity 
normalized by the speed of light, 7 = \/ ^Jl — (3'^ the bulk 
Lorentz factor, and R = ^/r'^ -\- a distance from the cen- 
tral object, In the present simple model, the mass-outflow 
rate M is assumed to be constant, while the wind velocity 
V is generally a function of R. 

We further assume that the heating and cooling may 
be ignored and the wind would expand adiabatically (cf. 
Park & Ostriker 2007). Then, the temperature Tg of the 
wind gas in the comoving frame is proportional to the 
density as Tq oc /Cq~^, where T is the ratio of specific heats 
and set to be 5/3 now. Hence, the comoving temperature 
varies as 

-4/3 

(2) 



To 



where Tc is the temperature at some reference radius Rc, 
that is fixed as Rc — r^. 

In this paper we assume a one-temperature plasma; i.e., 
the electron temperature To equals to the ion temperature 
Tq. If the flow is not adiabatic but there is energy ex- 
change between ions and electrons, the plasma generally 
becomes a two-temperature one. 

We suppose that the plasma is supplied from the hot 
gas accretion like ADAF in the equatorial plane including 
thermal and nonthermal electron populations. 

2.2. Thermal and Nonthermal Electron Populations 

Based upon a model by Ozel et al. (2000), we first fix 
the electron populations with a few parameters. The num- 
ber density iVth of electrons in the thermal population is 
assumed to be sufficiently larger than the number density 
iVpi of electrons in the nonthermal population: 
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and we set 
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Nth = Po/rup, 



(4) 



where mp is the proton mass. That is, it is assumed that 
the nonthermal electrons do not affect on the flow dy- 
namics and thermal properties, and a large fraction of the 
electrons in the flow are in a thermal distribution with 
temperature Te (= Tq). 

For the thermal electron population, we use the rela- 
tivistic Maxwell-Boltzmann distribution: 



?^th(7o) = A'th 



7o/3cexp(-7e/6'e) 



(5) 
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where 7e is the electron Lorentz factor, (3^ is the relativistic 
electron velocity, and 



(6) 



is the dimensionless electron temperature, and Ar2(l/^c) is 
the modified Bessel function of second order. The energy 
density of this Maxwell-Boltzmann distribution is derived 
as (Chandrasekhar 1939), 

tith a(6'c)A^tiimoC^6'c, (7) 

where a coefficient a{9c) is approximated as (Gammie & 
Popham 1998) 
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for the present purpose with sufficient accuracy. 

For the nonthermal electron population, on the other 
hand, we use a power-law distribution extending from 70 = 
1 to infinity: 



r^pl(7c)=A^pl(p-l)7e"^ 



(9) 



where the index p is a one of free parameters. The energy 
density of this power-law distribution becomes 



-fclTLcC^npid'-fc = 
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as long as p > 2. 

If we assume that a fraction 
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is constant in any radius, then the number density of the 
nonthermal electron is related to that of the thermal elec- 
tron by 



— ^ = rj- — -aiOcjOc 
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It should be noted that there is a mistype in equation (8) 
of Ozel et al. (2000). 

2.3. Synchrotron Emissivity 

Next, we summarize the synchrotron emissivity given 
by Ozel et al. (2000) (see also Rybicki & Lightman 1979; 
Mahadevan et al. 1996). 
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For the thermal electron population with the relativistic 
Maxwellian distribution, the synchrotron cmissivity jiy,th 
is approximated as 

Here, v is the frequency in the comoving frame, and the 
function M{x]\j) is well approximated by 



M{xm) 



4.0505 a 



1 



0.40 b 0.5316 c 



xexp(-1.8896a;^;'') 



1/4 
1/3n 



1/2 



where 
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The latter i^b is the nonrelativistic cyclotron frequency 
in a magnetic field of strength B, which is another free 
parameter. We set the coefficients a,b,c, which weakly 
depends on the temperature, asa = 6 = c= l for simplicity. 

For the nonthermal electron population with a power- 
law distribution, the cmissivity j^ ^i is given by (Rybicki 
& Lightman 1979) 

(l-p)/2 
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where 
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r(p/4 + 7/4) 

where F is the Gamma function. 

In contrast to the previous studies (e.g., Ozel et al. 
2000), we here examine the spectral modification by the 
relativistic effect. The observed frequency I'obs in the in- 
ertial frame is related to the emitted frequency v in the 
comoving frame by 



(19) 



where z is the redshift originated from the bulk motion of 
the optically thin black-hole wind and the gravitational 
redshift, although we ignore the light-bending in this pa- 
per (cf. Hutsemckers & Surdej 1990; Dorodnitsyn 2009 
for P Cyg profile of relativistic winds). Furthermore, the 
observed emissivities are also given as 



jL',th,obs 
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Finally, under the optically-thin assumption, the comov- 
ing and observed luminosities become, respectively. 
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In optically thin accretion flows thermal synchrotron 
and its Comptonization play an important role to the 
shape of the broad band spectrum (e.g., Kusunose & 
Takahara 1989; Narayan & Yi 1994; Mahadevan 1997; 
see also Park & Ostriker 2006). In particular, the 
Comptonized radiation is emitted in high energy spec- 
trum. However, the purpose of this study is not to fit the 
spectrum, but to demonstrate the relativistic effect in the 
outflow. Hence, we ignore the Comptonization. 

3. Results 

Using a wind model and approximated synchrotron 
emissivities, we calculate spectra from hot, optically thin 
black- hole winds for various parameters. 

The parameters of a wind model are the black hole mass 
M, the mass-outflow rate M, the wind (terminal) velocity 
/?, and the central temperature Tc, while the parameters 
of cmissivity are the fraction 77, the power-law index p, and 
the magnetic strength B. Of these, we fix M = lO^Af0 and 
Tc = 10^^ K, bearing in mind blazers and low luminosity 
active galaxies. The fraction of the nonthermal electron 
population is fixed as 77 = 0.04, using the recent results by 
Inoue et al. (2008). Hence, the rest parameters are /?, p, 
B, and to, where 



(24) 



is the mass-outflow rate normalized by the critical rate, 
Le being the Eddington luminosity of the central object. 

For wind velocity laws, we in turn examine two cases; 
constant and increasing cases. 

3.1. Constant Velocity 

We first consider the relativistic effect of the wind bulk 
motion on the observed spectrum, which is of most in- 
terest in this study. Spectra for typical parameters are 
shown in figure 1. 

In figure 1 the comoving and observed spectra are shown 
by thin dashed and thick solid curves, respectively, for 
several values of the wind velocity. The constant velocities 
are /? = 0.1,0.5,0.9 from top to bottom. Other parameters 
are fixed as m = 10^^^ p = 2.5, and B = 1 gauss. Here, 
the magnetic field strength B is assumed to be constant 
spatially. 

From figure 1, we first notice that the spectral intensity 
decreases as the wind velocity increases. This is because 
that the density of the observed region decreases as the 
wind velocity increases for the same m. 

Next, as the wind velocity increases, the observed spec- 
tral intensities become much more lower than the comov- 
ing ones, although the peak frequencies of the observed 
spectra shift toward the high frequency direction than 
those of the comoving ones. At first thought, the relativis- 
tic effect generally shifts the spectra toward a high energy 
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Fig. 1. Comoving spectra (thin dashed curves) and observed 
ones (thick soUd ones) for f3 = 0.1,0.5,0.9 from top to bot- 
tom. Other parameters are fixed as m = 10-3, p = 2.5, and 
B = I gauss. 



Fig. 2. Comoving spectra (thin dashed curves) and observed 
ones (thick solid ones) for p = 2.1,2.5,3 from top to bot- 
tom. Other parameters are fixed as /3 = 0.9, rh = lO"'^, and 
B = I gauss. 



regime. Indeed, in the optically thick case (Smnitomo 
et al. 2008), the observed spectra shift toward high fre- 
quency and high intensity directions, as the wind velocity 
increases. In the present case, the peak frequencies shift 
toward the high frequency direction, as expected. The 
spectral intensity, however, shift toward a low energy di- 
rection. 

This is understood as follows. In the present optically 
thin wind, we can observe the receding part of the flow and 
the part moving the perpendicular directions as well as the 
approaching part. When the flow velocity is sufficiently 
relativistic, the emitted intensities from the approaching 
part become blueshift, but those from the receding and 
perpendicular parts become redshift. As a result, the total 
observed spectra become lower than the emitted ones. 

Compared to the results by Ozel et al. (2000), the shape 
of the peak spectrum in figure 1 is somewhat flat. This is 
because we assume the uniform magnetic field. If we as- 
sume the equipartition, i.e., the gas pressure equals to the 
magnetic pressure, the magnetic field strength B is pro- 
portinal to r~^^^. Thus, the shape of the peak spectrum 
will be more sharpened. 

We further examine other parameter dependence of the 
present optically thin black-hole winds. 

In figure 2 the comoving and observed spectra are shown 
by thin dashed and thick solid curves, respectively, for 
several values of the index p. The values of p are p = 
2.1,2.5,3 from top to bottom. Other parameters are fixed 
as /3 = 0.9, m = 10""^, and B = 1 gauss. 

The steepness of the power-law part of the synchrotron 
emission is generally expected as L^, cx jy-^P-i)/^ (Rybicki 
& Lightman 1979). Since the Doppler effect does not 
change the power law steepness, the observed spectra shift 
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Fig. 3. Comoving spectra (thin dashed curves) and observed 
ones (thick solid ones) for B = 1,0.1,0.01 gauss from top to 
bottom. Other parameters are fixed as /3 = 0.9, rh = 10^^, 
and p = 2.5. 

toward the high frequency and low intensity directions in 
a self-similar manner, as the flow speed increases. 

In figure 3 the comoving and observed spectra are shown 
by thin dashed and thick solid curves, respectively, for 
several values of the magnetic strength B. The values of 
B are B ~ 1,0.1,0.01 gauss from top to bottom. Other 
parameters are fixed as /? = 0.9, rh = 10~^, and p = 2.5. 
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Fig. 4. Comoving spectra (thin dashed curves) and observed 
ones (thick soUd ones) for m = 10~^, 10~^, 10~* from top to 
bottom. Other parameters are fixed as f) = 0.9, p = 2.5, and 
-B = 1 gauss. 



Fig. 5. Comoving spectra (thin dashed curves) and observed 
ones (thick sohd ones) for /3oo = 0.1,0.5,0.9 from top to bot- 
tom. Other parameters are fixed as m = 10-^, p = 2.5, and 
B = 1 gauss. 



As the magnetic field becomes weak, the peak frequency 
of thermal parts shifts toward the low frequency direction 
and the intensity becomes lower. On the other hand, the 
relativistic bulk motion shifts the peak frequency toward 
the high frequency direction. 

In figure 4 the comoving and observed spectra are shown 
by thin dashed and thick solid curves, respectively, for 
several values of the mass-outflow rate m. The values of 
m are m = 10"2,10~3,10"^ from top to bottom. Other 
parameters are fixed as /3 = 0.9, p = 2.5, and B = 1 gauss. 

As is seen in figure 4, the mass-outflow rate does not 
affect the peak frequency, but only change the intensity. 
That is, as the mass-outflow rate decreases, the intensity 
simply decreases, since the gas density becomes low. As 
the flow velocity increases, the spectra shift toward high 
frequency and low intensity directions. 

3.2. Increasing Velocity 

In the realistic case, the hot gas must be accelerated 
from the very vicinity of the black hole to the outerside, 
where the wind speed would reach the constant terminal 
one. In this subsection, we thus consider such a realistic 
wind velocity law, and compare the previous case of the 
constant velocity law. 

As a velocity law, we adopt the monotonically increas- 
ing type, which is often used in the stellar wind: 

/3 = /3o + (/3oo-/3o)(^l-^) , (25) 

where the index a is set to be unity for simplicity. The 
initial radius and velocity are also set to be i?o = 2.0 Tg 
and /?o = 0.1, respectivey, while the terminal velocity (i^D 
is left as a free parameter. Therefore, the parameters are 



f3oc, P, B, and m. 

In figure 5 the comoving and observed spectra are shown 
by thin dashed and thick solid curves, respectively, for 
several values of the terminal speed (see also figure 1). 
The values of the terminal speed are /3oo =0.1,0.5,0.9 from 
top to bottom. Other parameters are fixed as m = 10^"^, 
p = 2.5, and B ~1 gauss. 

In the case of (3oo =0.1 the wind profile is constant and 
the spectra is the same as the case of = 0.1 in figure 1. 
On the other hand, in the cases of (3oo = 0.5 (and 0.9), the 
wind is accelerating from the law speed to high speed, and 
the results are different from those shown in figure 1. For 
example, the comoving spectra do not change so much, 
compared with the constant velocity case. That is, the 
relativistic effect is somewhat suppressed. This is because 
in the accelerating wind the velocity of the central part, 
where the synchrotron emissivity is large, is small. 

Moreover, the profiles of the thermal emission is dif- 
ferent from those in the constant velocity case. This is 
because the thermal emission strongly depends on the 
density A^th, which also strongly depends on the velocity 
profile. 

In figure 6 the comoving and observed spectra are shown 
by thin dashed and thick solid curves, respectively, for sev- 
eral values of the index p (see also figure 2). The values of 
p are p = 2.1,2.5,3 from top to bottom. Other parameters 
are fixed as f3oo = 0.9, m = 10"'^, and B = 1 gauss. 

As already stated, we notice that the spectral profile of 
the thermal part is rather different from the constant ve- 
locity case. Much more important is the total luminosity. 
Generally, in the realistic case with the increasing veloc- 
ity law, the luminosity decreases as the terminal speed 
increases since the density decreases. However, we found 



6 



N. Sumitomo et al. 



40- 




35- 



9 10 11 12 13 14 15 16 

log V 

Fig. 6. Comoving spectra (thin dashed curves) and observed 
ones (thick solid ones) for p = 2.1,2.5,3 from top to bot- 
tom. Other parameters are fixed as /3oo = 0.9, rh = 10""^, 
and B = I gauss. 

that the luminosity is rather higher than that of the con- 
stant velocity case. This is also because in the increasing 
velocity case the This is because in the accelerating wind 
the velocity of the central part is small. 

4. Concluding Remarks 

In this paper, we examined the observational appear- 
ance of hot, optically thin, rclativistic, spherically sym- 
metric black-hole winds from the observational point of 
view. We have calculated emitted and observed spectra 
of hybrid thermal-nonthermal synchrotron emissions, tak- 
ing into account the rclativistic effect in the high velocity 
regime. The spectral intensity generally decreases, as the 
wind velocity increases or the mass-outflow rate decreases. 

We could see deeper inside the wind, as the velocity 
increases. This is the relativistic limb-darkening effect. 
This nature does not depend on the observer's direction. 
In addition, the luminosity in the observer's frame is re- 
markably enhanced by relativistic beaming effects along 
the observer's direction. These two effects mainly work 
as the luminosity enhancement of the rclativistic outflow. 
We suggest that if the observed luminosity is used for the 
evaluation of the black hole mass, then the derived black 
hole mass will be overestimated. 

We found that the intensity of the observed spectra de- 
creases compared with that of the emitted ones, although 
the peak frequency shifts towards the high frequency di- 
rection. For example, in the cases of /3 = 0.9 or /3oo — 0.9 
the observed intensity becomes smaller than the emitted 
one more than about one order. This is because in the 
optically thin outflows we can see the far side flows that 
go away from the observer. 



Furthermore, we found that the intensity of the con- 
stant velocity case becomes smaller than that of the ac- 
celerating case. This is because in the accelerated velocity 
case the central part has a low velocity. We thus carefully 
consider optically thin relativistic flows around a black 
hole such as Sgr A*. 

In the present study we did not solve the radiative trans- 
fer equation, while Ozel et al. (2000) solved the radiative 
transfer to obtain the hybrid synchrotron spectra. We fo- 
cus our attention on the relativistic effect of hot, optically 
thin black-hole winds, and the relativistic effect is gener- 
ally important in the high frequency regime. In order to 
solve the spectra in the low frequency regime correctly, we 
should solve the radiative transfer equation as Ozel et al. 
(2000) have done. 

Moreover, we ignored the general relativistic bending 
of light, which becomes important at the very center of 
winds (Hutsemekers & Surdej 1990; Dorodnitsyn 2009 for 
P Cyg profile of relativistic winds). In the realistic stua- 
tion winds would be accelerated at the very center. We 
should include these detailed behavior at the very center 
of winds in future. 

Finally, we ignored the Comptonization in this paper 
(e.g., Kusunose & Takahara 1989; Narayan & Yi 1994; 
Mahadevan 1997; Park & Ostriker 2006). In the hot 
plasma, however, the emitted bremsstrahlung and syn- 
chrotron photons are generally upscattered by inverse 
Compton scattering. For example, the Compton-heated 
outflow was examined in details in Park and Ostriker 
(2006). Such a Comptonization works for the gas as 
cooling, and modifies the spectra (Rybicki and Lightman 
1979). Although in order to demonstrate the relativistic 
effect we ignored the Comptonization in this paper, we 
should include the Comptonization to give more realistic 
spectra from the hot black- hole wind. 

References 

Abramowicz, M. A., Novikov, I. D., & Pacyriski B. 1991, ApJ, 
369, 175 

Begelman, M. C. & Chiueh, T. 1988, ApJ, 332, 872 
Bisnovatyi-Kogan, G. S., & Lovelace, R. V. E. 1997, ApJ, 486, 
L43 

Colpi, M., Maraschi, L., & Treves, A. 1984, ApJ, 280, 319 
Dorodnitsyn, A.V. 2009, MNRAS, 393, 1433 
Fukue, J., & Sumitomo, N. 2009, PASJ, 61, in press 
Gammie, C.F., & Popham, R. 1998, ApJ, 498, 313 
Ghisellini, G., Haardt, F., & Fabian, A. C. 1993, MNRAS, 263, 
L9 

Ghisellini, G., Haardt, F., & Svensson, R. 1998, MNRAS, 297, 
348 

Gruzinov, A., & Quataert, E. 1999, ApJ, 520, 849 
Hutsemekers, D., & Surdej, J. 1990, ApJ, 361, 367 
Inoue, Y., Totani, T., & Ueda, Y., 2008, ApJ, 672, L5 
Kato, S., Fukue, J., & Mineshige, S. 2008, Black-Hole 

Accretion Disks (Kyoto: Kyoto University Press) 
Kusunose, M., & Takahara, F. 1989, PASJ, 41, 263 
Li, H., Kusunose, M., & Liang, E. P. 1996, ApJ, 460, L29 
Mahadevan, R., Narayan, R., & Yi, I. 1996, ApJ, 465, 327 
Mahadevan, R. 1997, ApJ, 477, 585 
Medvedev, M. 2000, ApJ, 541, 811 



Relativistic Spectra of Hot Black-Hole Winds 



7 



Narayan, R., & Yi, I. 1995, ApJ, 452, 710 

Narayan, R., MaClintock, J.E., Yi, I. 1996, ApJ, 457, 821 

Ozel, F., Psaltis, D. & Narayan, R., 2000, ApJ, 541, 234 

Park, M-G., & Ostriker, J. P. 2007, ApJ, 655, 88 

Quataert, E., & Gruzinov, A. 1999, ApJ, 520, 248 

Quataert, E., & Narayan, R. 1999, ApJ, 520, 298 

Rybicki, G.B., & Lightman, A. P. 1979, Radiative Processes in 

Astrophysics (New York: John Wiley & Sons) 
Sumitomo, N., Nishiyama, S., Akizuki, C., Watarai, K., & 

Fukue, J. 2008, PASJ, 59, 1043 
Zdziarski, A. A., Lightman, A. P., & Maciolek-Niedzwiecki, A. 

1993, ApJ, 414, L93 



